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We analyze structural and conformational properties in a simulated bead-spring model of a non- 
entangled, supercooled polymer melt. We explore the statics of the model via various structure 
factors, involving not only the monomers, but also the center of mass (CM). We find that the 
conformation of the chains and the CM-CM structure factor, which is well described by a recently 
proposed approximation [Krakoviack et al, Europhys. Lett. 58, 53 (2002)], remain essentially un- 
changed on cooling toward the critical glass transition temperature Tc of mode-coupling theory. 
Spatial correlations between monomers on difi'erent chains, however, depend on temperature, albeit 
smoothly. This implies that the glassy behavior of our model cannot result from static intra-chain or 
CM-CM correlations. It must be related to inter-chain correlations at the monomer level. Addition- 
ally, we study the dependence of inter-chain correlation functions on the position of the monomer 
along the chain backbone. We find that this site-dependence can be well accounted for by a theory 
based on the polymer reference interaction site model (PRISM). We also analyze triple correlations 
by means of the three-monomer structure factors for the melt and for the chains. These structure 
factors are compared with the convolution approximation that factorizes them into a product of two- 
monomer structure factors. For the chains this factorization works very well, indicating that chain 
connectivity does not introduce special triple correlations in our model. For the melt deviations are 
more pronounced, particularly at wave vectors close to the maximum of the static structure factor. 

PACS numbers: 61.25.Hq Macromolecular and polymer solutions; polymer melts, 61.20.Ja Computer simu- 
lation of liquid structure 



I. INTRODUCTION 



The microscopic understanding of the glass transition 
is a challenging problem in contemporary condensed mat- 
ter physics 00,0. During the past decade the re- 
search in this field was stronjgly influenced by the mode- 
coupHng theory (MCT) H U This theory suggests that 
a non-linear coupling between density fluctuations drives 
the slowing down of the structural relaxation when a li- 
quid approaches its glass transition. MCT predicts that 
there is a critical temperature Tc, experimentally found 
to be above the calorimetric glass transition tempera- 
ture Tg, where the dynamics qualitatively changes. For 
T > Tc the relaxation of the glass former is determined by 
the cooperative motion of the particles comprised in the 
nearest- neighbor shells ( "cage-effect" ) . In "ideal MCT" , 
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the simplest version of the theory, the mutual blocking of 
the particles in the cages leads to a complete structural 
arrest at T^. This complete freezing is not observed ex- 
perimentally, possibly with the exception of polydisperse 
hard-sphere- like colloidal suspensions 0. When cooling 
the glass former toward Tg the structural relaxation time 
continuously increases, instead of diverging at Tc. Thus, 
alternative relaxation mechanisms must exist besides the 
cage effect and eventually become dominant for T < Tc. 
Within MCT the microscopic origin of these processes is, 
however, not well understood. 

Despite this limitation of its range of validity the ideal 
MCT has been tested in numerous experiments and 
computer simulations 0,0,01 ■ Broadly speaking, the the- 
oretical predictions were found to provide an adequate 
description of the relaxation dynamics above Tc. This 
success has stimulated extensions of the theory, originally 
developed for simple liquids, to molecules with orienta- 
tional degrees of freedom 0, ^ and recently also to 
polymers 
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A distinguishing feature of the theory consists in es- 
tabhshing a quantitative Hnk between the structure of a 
glass former and its dynamics. This link may be exploited 
to predict the relaxation behavior provided the relevant 
static properties are available. These properties involve 
the static structure factor and related quantities which 
must be determined with high precision over a large range 
of wave vectors. Presumably due to this prerequisite a 
quantitative comparison of the predicted and measured 
dynamics has been attempted only for a few systems in 
the past, such as hard-sphere- like colloidal particles [l^ . 
soft-s phe re [T^ . hard-sphere Il4l or Lennard- Jones mix- 
tures [13 , diatomic molecule s 11(1. a nd models for ortho- 
terphenyl [l^ and Si02 |l9l l2(]|| . These studies sug- 
gest that for T > Tc MCT is a promising approach to a 
quantitative description of the structural relaxation for a 
large class of liquids comprising fragile and strong glass 
formers. 

These findings motivate our present work. An exten- 
sion of MCT to non-entangled polymer melts |^ |^ 
opens the possibility to attempt a quantitative compar- 
ison also for a polymeric glass former. Here, we present 
the first step toward such a comparison for a simu- 
lated bead-spring model of a supercooled polymer melt 
[22. l2a. KM. I25I l2g | . We discuss various static structure 
factors, paying particular attention to the dependence of 
the structure on the position of a monomer along the 
polymer backbone and to correlation functions involving 
the center of the mass of the chains. This information 
may be used to develop a tractable theory. We will re- 
port on that and on the comparison with the simula- 
tions in a forthcoming article [2ll . A key aspect of this 
theory is that the short-range order of the monomers, 
as measured by the main peak of the collective static 
structure factor, strengthens with decreasing tempera- 
ture. The strengthening of the local packing provides 
the dominant mechanism causing structural arrest and 
glassy dynamics. This mechanism, termed "cage effect" 
in simple liquids is also at the core of our theory for 

polymer melts. In this article, we will thus pay special 
attention to structural correlations around the average 
spacing between monomers, and to the question of how 
chain connectivity affects them. 



being the minimum of the non-truncated potential. 

In addition to C/lj, successive monomers along 
polymer backbone interact via a FENE potential |21 
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with Rq — 1.5 and k = 30. The superposition of the 
LJ- and FENE potentials leads to a steep effective bond 
potential with a sharp minimum at rb = 0.9606. 

This choice of parameters has two important conse- 
quences. First, it prevents bonds from crossing each 
other. This imposes topological constraints j2^] which 
ultimately lead to reptation-like dynamics in the limit 
of long chains [2?! l2^. Second, the bond potential lo- 
cally distorts the regular arrangement of the monomers 
because it favors the inter-monomer distance rb which is 
incompatible with r„iin- 

When cooling the melt from high T the incompatibility 
of Tb an d rrn jn iinpedes crystallization, but does not pre- 
clude it [33, IsJ, |32| ■ For the melt to remain amorphous 
the chains should also be flexible. This was pointed out in 
simulations of a semi-flexible bead-spring model in which 
large bond angles are energetically favored by a bending 
potential The interplay of chain stiffness and 

excluded volume interactions suffices to induce crystal- 
lization from the melt. Contrary to that, the chains of 
our model are flexible. In the temperature range stud- 
ied, the end-to-end distance {R^ ~ 12.3) and the radius 
of gyration (i?g ~ 2.09) are almost constant, and the col- 
lective static structure factor of the melt is typical of an 
amorphous material pMl25|. 

We analyze time series of isobaric simulations at the 
pressure p = \ [s^l • The polymer melt contains n 
monodisperse chains of length = 10 in the volume V . 
Depending on temperature (Nose-Hoover thermostat) n 
ranges between n = 100 and n — 120. This corresponds 
to the following chain (p) and monomer densities (pm) 



0.091 < p = — < 0.104 , 

nN 

0.91 < Pm = -^^ < 1.04 



(3) 



II. MODEL AND SIMULATION TECHNIQUE 

We study a bead-spring model of linear polymer chains 
HillS' AH monomers interact via a truncated and shifted 
Lennard-Jones (LJ) potential 

, Ae[{a/rf^~{a/rf]+C, r < 2r^^ , 

^ 1 0, r>2rmi„. 

(1) 

In the sequel, we will use LJ units (e = 1, ct = 1; further- 
more, Boltzmann's constant — \ and the monomer 
mass m = 1). The constant C — 127/4096 is chosen so 
that JTlj vanishes continuously at r = 2rniin, rmin = 2^/^ 



III. THEORETICAL BACKGROUND 

A. Basic Notations 

Let r° denote the position of the ath monomer in chain 
i and the position of the center of mass (CM) of chain 

1 ^ 

^'^=nY.< (^-1, ■•■,«)• (4) 

The knowledge of and R, allows us to define various 
density fluctuations for the wave vector q in reciprocal 
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space: the density fluctuations of monomer a, 

n 

Pa(q) -5]cxp[iq.<] (a=l,...,7V), (5) 

the density fluctuations of a tagged chain i obtained by 
summing over all monomers of the chain, 

N 

pP(q)=^exp[iq.ra , (6) 

a=l 

the density fluctuations created by all monomers of the 
melt, 

n N 

Ptot(q) =^^exp[iqT,"] , (7) 

i—l a—1 

and the polymer-density fluctuations related to the CM's 
of all chains 

n 

1=1 

Density-density correlation functions are well established 
means to describe the structure (and the dynamics) of a 
liquid 33] . For a polymer melt we can derive various such 
two-point correlation functions from Eqs. ||SHHI)- They 
are introduced in the following section. 

B. Static Structure Factors 

The density-density correlations at the monomer level 
may be characterized by the monomer-monomer (or site- 
site) static structure factors 

Sabiq) = - (/Oa(q)*Pb(q)) • (9) 
n 

Here, (•) denotes the canonical average over all configura- 
tions of the melt. Since the melt is spatially homogeneous 
and isotropic, the structure factors depend only on the 
modulus of the wave vector, |q| = q. 

We can split Eq. ((UJ) into an intra-chain and an inter- 
chain part: 

Sab{q) = Wab{q) + phab{q) ■ (10) 

The intra-chain contribution is given by 
1 " 

^-''('?) ^ - (Ee^p{-iq • - ^']}) ' (11) 

1=1 

and the inter-chain contribution by 

n 

phabiq) = -(^exp{-iq. [r,?-r5]}) . (12) 




FIG. 1: Schematic illustration of the correlation functions 
defined in Sect. 1111^ r" is the position of the ath monomer of 
chain i and Ri the position of the chain's center of mass. Wcd 
and 'Wa,c [Eqs. lllllSU l denote intra-chain structure factors, 
hcd and ha.c [Eqs. p2ll9|l ] inter-chain structure factors. Sc 
[Eq. 1161 1 is the structure factor of the centers of mass. 

These contributions reveal static correlations between 
monomers belonging to the same chain or to different 
chains (cf. Fig.P). 

When averaging over all monomer pairs (a, 6) we ob- 
tain the collective static structure factor of the melt, 

1 1 ^ 

Siq) = ^ (Ptot(q)>tot(q)) = ]^ $Z 

a,h=l 

= W{q) + PnMq) : (13) 
where pm — Np and 

w{q) = - (pr(q) vr(q)) = ]^ E ^^^^(9) (14) 

a, 6=1 

denotes the static structure factor of a chain and 

^^(9) = ]^ E ^-f(9) (15) 

a, 6=1 

is the Fourier transform of the site-averaged inter- 
molecular pair-correlation function l33l|. 

Usually, the averaged quantities S{q) and w{q) are used 
to characterize the structure of a polymer melt. Contrary 
to that, we focus here on the monomer-resolved genera- 
lizations Sabiq) ^-iid Wabiq)- The aim of our study is to 
understand to what extent specific monomer- monomer 
correlations deviate from the average behavior. Since 
the structure factors Sabiq) "Wabiq) are important in- 
put quantities for the mode-coupling approach to glassy 
polymer dynamics 11], the comparison of Sabiq) and 
Wabiq) with their monomer- averaged counterparts can 
suggest suitable approximations and thus help develop- 
ing a tractable theory . 

In addition to density fluctuations of the monomers 
the spatial arrangement of the CM's and the coupling 
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between the CM and the monomers can be analyzed. We 
define the CM-CM structure factor (see Fig. [TJ 



Sciq) = -(pc(q)*pc(q)) 
1 " 

= - {Yl exp{-iq- [R, - R^]} 



= 1 + phc{q) , 



(16) 



which we spht, in analogy to Eq. into self- ("1") 

and distinct (he) parts. Formally, Eq. Hlt)|) is identical 
to that of simple hquids [ssf . 

Similarly, the coupling between a monomer and the 
CM's leads to monomer-polymer structure factors: 



1 



Safiiq) - -(Pa(q)>c(q)) 



1 " 

-(5:exp{-iq.[r^-R,]} 



(17) 



Sa,ciq) is Fourier transform of the (averaged) prob- 
ability of finding a site a at a distance r from the center 
of a chain. Following Eqs. H1UH12|I we separate again the 
intra-chain contribution. 



1 " 

Wa,c (?) = -( XI exp{ -iq • [r^ - R;] } 



(18) 



from the contribution involving different chains 
1 " 

pha,c{q) ^ -{T.^'M-^^- i^"} -^j]}) ■ (19) 



The correlations measured by Eqs. H18I19|) are illustrated 
in Fig.^ While the intra-chain structure factor Wa,c(<z) 
can also be determined for a single polymer at infinite di- 
lution '34'! , the inter-chain ha,c{q) describes how a-sites 
arrange around the CM of another polymer. When sum- 
ming over all monomers and defining 



N ^ N 



(20) 



we obtain from Eqs. ((TTHT^ 

N 

Sm,c{q) = X! '^a:C(9) = Wni,c{q) + PmV,c(g) • (21) 



These functions describe the averaged packing of sites 
around the center of mass of the same (i«rn,c) a-nd of a 
different (/im,c) polymer, respectively. 



C. Three-Particle Structure Factors 

S'(q) measures the spatial distribution of two 
monomers. It depends on one wave vector q. Its gen- 
eralization to a triple correlation function, the three- 
monomer structure factor 6*3 (q,k), describes the (aver- 
aged) packing of a third monomer which results from 
fixing the position of two monomers. 53(q, k) depends 
on two wave vectors and is defined by 



'S'3(q,k) = ^(ptot(-q)ptot(k)/9tot(q - k)) = 



1 



n N 
— l a,b,c—l 



exp{i[-q • + k • + (q - k) • r^} 



(22) 



where the vector q — k = p is the third side of a triangle 
forming an angle f between q and k. The angle ip is 
given by 



cos If 



q 



2qk 



(23) 



5*3 (q,k) is the Fourier transform of the three particle 
distribution function, which gives the average density of 
(other) segments at a space point if the positions of two 
segments are fixed. It is related to a direct correlation 
function C3 (q, k) by the triplet Ornstein-Zernike equation 

M 

53(q, k) = S{q)S{k)S{\q - k|) [l + pldq, k)] . (24) 

Often, when considering 6*3 for wavevectors correspond- 
ing to the average particle distance, it is assumed that 



r 



there are no three-body correlations which are not con- 
tained in the product of two-particle correlation func- 
tions. This so-called convolution approximation 35, 3^] 
-note that it differs from the Kirkwood superposition ap- 
proximation which becomes valid for large distances |33l|- 
implies that the triplet direct correlation function van- 
ishes, C3(q, k) = 0. Here, we want to test this approxi- 
mation for 53 (q, k) and for the polymer three-monomer 
structure factor W3(q, k). ?i'3(q, k) is defined analogously 
to S3 with ptot replaced by in Eq. H22I) . that is, only 
the term i = j = I oi the sum in Eq. H22|l is taken into 
account. 

For a homogeneous and isotropic system depends 
on the moduli of the three wave vectors only, 53 (q, k) = 
53(9, k,p = |q — k|). To determine S'3((7, k,p) we utilized 
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a method closely related to the one proposed in Ref. [20j . 
The triple of moduli {q,k,p) satisfying Eq. is dis- 
cretized in bins of width Ag = A A; = Ap = 0.2. In each 
bin 100 vector tuples {(q, k) | |q| = q, |k| — k, |p| — p)} 
are chosen at random for each configuration and S{q, k,p) 
is calculated as the average over this set of vectors and all 
configurations. Data were accumulated over 1155 config- 
urations at T = 0.47 IS. 



IV. RESULTS ON TWO-POINT CORRELATION 
FUNCTIONS 

A. Static Properties at the Monomer Level: 
Site-averaged Quantities 

We discuss the static structure factor of the melt and 
the corresponding self- and distinct parts [Eqs. (|10m5() ]. 
For T < 1 and q < 20, results for SYg) and w{q) have 
been presented previously 0,123, |2l,|2g. Here, we extend 
the analyses up to g = 50 and considerably improve the 
statistics (averages over more than 1000 configurations). 
This effort was necessary to use the static quantities in 
MCT calculations which require a large g-range and good 
statistics. 



1. Static Structure Factor of the Melt 

Figure 13 shows S{q) for temperatures above the criti- 
cal temperature of MCT (Tc ~ 0.45). In this T-interval 
the structure of the melt is typical of a disordered, dense 
system. Due to the weak compressibility of the melt S{q) 
is small in the limit q 0. As q increases, S{q) increases 
toward a maximum which occurs around gmax — 7.15 
in our model. This g-value corresponds to the length 
scale of the monomer diameter (= 1). Thus, the domi- 
nant contribution to S'(gmax) comes from the amorphous 
packing in the nearest neighbor shell around a monomer. 
On cooling, Fig. |21 indicates that no long-range struc- 
tural correlations develop in the melt. Only its density 
increases and the packing becomes tighter. Two features 
refiect these changes: the amplitude of the peak S'(gmax) 
grows and its position gmax shifts to larger values. Close 
to gmax, the dependence of S{q) on T is most pronounced. 

The increase of S'(gmax) with decreasing T may be re- 
lated to the (empirical) Hansen- Verlet freezing criterion 
|39l |. This criterion states that a liquid will condense to 
a solid phase when S'(gmax) exceeds a threshold value. 
For crystalhzation this threshold is ^(gmax) ~ 2.85 |40l |. 
For the glass transition of hard spheres within the ideal 
MCT one finds S'(gmax) ~ 3.54 "s^. The latter value 
agrees very well with our simulation result for S'(gmax) 
at T = 0.46 (see inset of Fig.|21), indicating that the melt 
is close to the T^oi the ideal MCT. Indeed, our previous 
studies mmillli,!!! suggest ~ 0.45. 

This interpretation of the structure factor only involves 
packing arguments which could also be put forward for 
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FIG. 2: Collective static structure factor S{q) of the melt 
versus the modulus of the wave vector q [Eq. CSJ]. The 
temperatures shown are: T = 1 (dashed line), 0.7, 0.65, 0.6, 
0.55, 0.52, 0.5, 0.48, 0.47, and 0.46 (solid line). Inset: First 
maximum of S{q), S{qmax), versus T (gmax ~ 7.15). The 
dashed horizontal line indicates the Hansen- Verlet criterion 
for the glass transition of hard spheres within the ideal MCT 
(5(gn.ax) ~ 3.54) p. 



atomic liquids. To obtain a better insight into the role of 
chain connectivity we split S{q) into intra-chain (self) 
and inter-chain (distinct) contributions [Eqs. (|10H15|) ]. 
We discuss these contributions in the next two sections. 



2. Intra-chain Structure 

Figure 131 shows the intra-chain contribution to 5'(g), 
the structure factor w{q) of a polymer, at T = 0.46 and 
T = 1. Both temperatures yield almost identical results 
for all g. So, data at intermediate T are not included in 
the figure. 

The independence of T may be rationalized in the fol- 
lowing way. If g ^ oo, local rapid variations will deter- 
mine the behavior of w(g). In this limit, we can approx- 
imate Eq. l(TI|l by 



w{q) 



l + —(N 



l)waa+i{q) (9 large) 



(25) 



where we also assumed that the nearest-neighbor con- 
tribution Waa+i{q) is independent of a. This assump- 
tion is well justified for our model, as a glance at Fig. El 
shows. Waa+iiq) can be evaluated by exploiting that 



Thus, 



rb due the stiff bond potential of our model. 



Wa a+liq) 



sin(g[r° — r°+^ [) \ sin(grb) 



g|r" 



qrh 



To lowest order, we therefore expect 



2 , , sin(grb) , , 

w{q) ~ I + ^AN ~ 1) ^Uii (g large) 
N qrh 



(26) 



(27) 



FIG. 3: Polymer static structure factor io(q) [Eq. II14II 1 at the 
lowest and the highest temperature, T — 0.46 (solid line) and 
T — 1.0 (dashed line), respectively. The simulation data are 
compared to the Debye formula (Eq. 12811 . dotted line) and to 
the large-q approximation Eq. 12711 (rb = 0.9609, dash-dotted 
line). The inset magnifies of the small-g behavior of w{q). 



FIG. 4: Distinct contribution to the static structure factor 
prnh{q) [Eq. ^] at T = 0.46, 0.65 and 1.0. Inset: Compar- 
ison of prah{q) with S{q) — 1 and — uid(ij) at T = 0.46 (see 
text for further details). 



For small g-values polymer physics suggests that the 
Debye function wj) (q) [23, provides a reasonable de- 
scription of w{q), 

wu{q)^Nfo{q^Rl), /d(x) = A[e-+a;-l] . 

(28) 

The Debye theory assumes a Gaussian distribution for all 
distances along the backbone of the chain. This assump- 
tion correctly reflects the random-walk-like structure of 
(long) polymers on intermediate and large length scales 
in the melt. However, at large q where the precise form 
of the interaction potential matters, it cannot apply. 

The approximations H27(l and (|28|l are compared to 
the simulation data in Fig. |2| For small q, the scat- 
tering is determined by the size of the polymer only. In 
this limit, the Debye theory reproduces the exact result 
'"^d('7) = .^^(1 — q^Rg/3). So, it has to coincide with the 
measured w{q) if q < 1/Rg — 0.69 (see inset of Fig. OJ. 
Furthermore, the theory and the simulation should also 
agree for l/i?g ^ q <C l/^b. This g-range probes the 
random-walk-likc internal structure of a polymer, where 
both w{q) and w-£i{q) scale as ~ q~^- However, since 
our chains are short, the length scales rb and _Rg are not 
sufficiently separated for this behavior to be observed. 
Instead, w{q) crosses over to regular oscillations which 
are compatible with Eq. H27|) for > 8. 

Figure 131 shows that the superposition of Eqs. 1)27(1 and 
(I28|l approximately describes the simulation data. Thus, 
in our model the main features of w(g) are determined 
by two length scales, rb and i?g. Since i?g depends only 
weakly on T and rb is independent of T, the almost per- 
fect agreement of the structure factors for T — 0.46 and 
T = 1 can be understood. 



3. Inter- chain Structure 



If the intra-chain contribution to S{q) is independent of 
T, the temperature dependence of S{q) must result from 
inter-chain correlations. Figure 0] supports this expecta- 
tion. The distinct part p^h{q) [Eq. (|13|l ] exhibits liquid- 
like oscillations whose extrema become more pronounced 
on cooling. This trend is especially visible around gmax 
(~ 7.15), suggesting that the glassy behavior of our 
model is driven by h{q) only. That is, by nearest neigh- 
bors which are not bonded to one another. This finding 
is not unreasonable. As the distance between successive 
monomers along the backbone is almost fixed, only non- 
bonded neighbors can pack more tightly and reinforce 
the cage with decreasing T. 

We support this interpretation by the following argu- 
ment: If our system was a simple liquid, we would have 
w{q) — 1, and p-^h{q) — S{q) — 1 would be exact. The 
inset of Fig. 21 shows that this simple-liquid-like appro- 
ximation represents a good description for 9^6. For 
smaller 5, deviations are found. Here, S{q) becomes van- 
ishingly small [S{q) ~ 10~^, see Fig. and polymer- 
specific effects, i.e., the correlation hole determine 
the distinct part: p,^h(q) « —w{q). The correlation hole 
implies that the probability of finding monomers of other 
polymers inside the volume Vc occupied by a chain is 
decreased. This effect arises because the probability of 
finding the monomers of a chain in its own volume Vc 
is enhanced, and intra- and inter-moleculer correlations 
have to compensate each other to render the melt incom- 
pressible [iali^. 
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B. Static Properties at the Monomer Level: 
Site-resolved Quantities 

1. Direct Correlation Function 

For simple liquids the direct correlation function is 
usually introduced through the Ornstein-Zernike equa- 
tion j^]. For molecular liquids Chandler and Ander- 
sen jZJ extended the Ornstein-Zernike approach to re- 
flect the contribution from the intra-molecular correla- 
tions WabiQ)- The resulting generalized site-site Ornstein- 
Zernike equation -also referred to as "Reference Interac- 
tion Site Model" (RISM)- is given by jH 

N 

habiq) ^ ^ Wax{q) Cxy{q) [wyb{q) + phyb{q)] . (29) 



Here, Cab{q) is the direct correlation function between the 
sites a and h. Inserting Eq. (|29|l into Eq. (|10|) we obtain 



pcab{q) = [wj{q) - Sj{q)\ 



(30) 



where X^^{q) denotes the (a, 6)-element of the inverse of 
the matrix X((7). 

The difficulty in dealing with site-site correlation func- 
tions arises from the dependence on the indices (a, fe). 
This leads to O(iV^) coupled equations which cannot 
be handled for large N . However, one can argue that, 
for long chains, end effects on inter-polymer correlations 
should be small, suggesting to treat all sites of a ho- 
mopolymer equivalently. (This simplification becomes 
exact for ring homopolymers.) This equivalent- site ap- 
proximation is usually invoked for Cab{q), i-e., 

Cab{q) = c{q) (equivalent-site approximation) . (31) 

Equation represents the principal idea of the poly- 
mer RISM (PRISM) theory developed by Schweizer and 
coworkers 43]. 

Inserting the assumption (|31|l into Eq. H29|) we obtain 
from Eq. H15|) the so-called PRISM equation for the site- 
averaged pair-correlation function 



h{q) = w{q) c{q) [w(g) -f p^Hq)] 



(32) 




Equations H13I32I) provide an expression for c{q) in terms 
of w{q) and S{q): 



FIG. 5: Examination of the equivalent-site approximation 
Eq. II31II at T = 0.47. The solid lines in both panels denote 
c(q) determined from the simulation results for S{q) and 'w{q) 
via Eq. 13311 . The dashed lines represent cii(g) (upper panel) 
and Ci5(g) (lower panel), the dotted lines 022(9), 233(5), 244(9), 
C55(g) (upper panel) and 025(9), £35(9), 045(9), Cssiq) (lower 
panel). The site-site direct correlation functions Cabiq) are 
calculated from the simulation results for Sab{q) and Wab{q) 
via Eq. 13UII . The insets magnify the region close to 9max. 



for the dynamics of the melt can be derived by assum- 
ing Eq. H31|) without introducing a large error [23 • In 
the following, we want to use Eq. (|31|l to interpret the 
inter-molecular site-site-correlations of our model. 



2. Intra-chain Structure 



Pmciq) 



w{q) S{q) 



(33) 



Figure |S1 examines the validity of the equivalent-site 
approximation by comparing Eq. (|30|l and Eq. I|33|) . Ap- 
parently, the approximation is well satisfied, except for 
functions involving the chain ends (a = 1 or a = N). 
Here, we find slight deviations close to gmax and more 
pronounced ones for g ^ 5. From the point of view of 
MCT, the important wave-vector regime is around gmax- 
Thus, Fig.|Slsuggests that, for our model, MCT equations 



If w{q) does not depend on temperature, we may ex- 
pect the same feature for all components Wab{q)- Indeed, 
we find that Wab{q) is (almost) independent of T. Thus, 
we concentrate on one temperature in the following. 

Figure ini depicts Wab{q) at T = 0.47 for various mono- 
mer pairs (a, 6). We find that the explicit dependence of 
Wab{q) on the site indices is negligible for all wave vectors. 
To a very good approximation, Wab(q) only depends on 
the distance |a — &|, a feature also found for a model of 
freely-jointed hard-sphere chains ^4^ and for the Gaus- 
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FIG. 6: Static structure factor Wabiq) of the monomer pair 
(a, fe) at T = 0.47. a and b are monomers of the same chain. 
Note that Wabiq) depends on |ffl — fe| only. The simulation data 
for \a — b\ = 1 {a = l,b = 2 and a = 5, 6 = 6) are compared 
with Eq. 1261 (circles). For separations |a — &| = 1,...,5 the 
Gaussian approximation Eq. 1341 is also shown (dotted lines) . 



sian approximation 



exp 



2| 

q \a- 



(34) 



To compare Eq. H34I) with the simulation data we iden- 
tify the statistical segment length a with cr'^ — R^/N — 
6Rg/N |2^. This assumption is valid for (long) chains 
in the melt. Figure shows that Eq. H34|) only provides 
a good approximation for small q. (The same result was 
also found when comparing w{q) and wu{q)-) With in- 
creasing q, Wab[q) decays faster than w^j^{q)^ if |a— 6| > 1, 
and becomes negative before approaching from below. 
This undershoot shifts to larger q and increases in ampli- 
tude, as |a — 6| decreases toward 1. Adjacent monomers 
along the backbone of the chain (|a — 6| = 1) exhibit long- 
range oscillations which are well described by Eq. 



3. Inter-chain Structure 

Figure [7| shows the site-resolved pair correlations 
hab{q) [Eqs. CUa)] at T 0.47. For q > 15, habiq) 
is independent of the monomer index and coincides with 
the site-averaged h(q). Contrary to that, we find a depen- 
dence on (a, b) around Q'max and particularly for q < 5.5. 
While pairs resulting from inner sites (1 < a, 6 < TV), ex- 
cept h22{q), are still very close to each other, differences 
occur for correlations comprising an end monomer. For 
g ~ 4, habiq) exhibits a shoulder. If an end monomer is 
involved, the amplitude of the shoulder decreases and is 
smallest for the correlation of two chain ends (hn). 

Qualitatively, the site-dependence around gmax may be 
explained by the following argument. As the analysis in 
Sect. lIV C 21 will show, a middle segment is buried deeply 
in the polymer coil and is closer to the CM of its chain 



FIG. 7: Inter-chain static structure factor, phab{q), at T = 
0.47 for different pairs (a, 6) [Eq. 1121 1. phab{q) depends on 
{a,b) at gmax — 7.15 and for q < 5.5. The correlation of two 
chain ends {a — b = 1) behaves differently in comparison to 
all other curves. The average over all monomers, ph{q), is 
also included. 



than is an end segment. Thus, it is not surprising that 
end monomers have stronger local inter-molecular corre- 
lations, viz. that phii{q) deviates more from zero around 
(Zmax than does phaaisi) for middle segments. 

A quantitative explanation of the site-dependence 
of habiq) ™ay be obtained by PRISM theory. From 
Eqs. (Iinj) and ^ it follows that 



Sab{q) = [{i-pHq)c{q)}-'Hq)]ab' 

habiq) = [w(g)c((7){I-pw(g)c((7)}"^w(g)]^j^ 



(35) 



where I denotes the unit matrix. Thus, even with the 
assumption Cabiq) — c(q), a site-dependence of Sabiq) 
and habiq) I'esults from chain connectivity due to the 
matrix structure of Wabiq)- E.g., for habiq) oi^^ finds 
[from Eqs. ^^M^] 



habiq) 



hjq) 
wiq)^ 



N 



^Waxiq) 



r N 
^Wbyiq) 

y=l 



(36) 



Figure IHl compares theory and simulation for some repre- 
sentative pair-correlation functions habiq)- We find that 
habiq) is vfG^^ described by Eq. H36|l. This explains why 
the correlation hole in habiq) for end segments is slightly 
more narrow than for middle segments. Furthermore, ac- 
cording to Eq. (|36|l hab iq) should exhibit the symmetry: 
ha.biq) — ha.N-b+iiq)- This Symmetry is tested in Fig.|S| 
for a — 2, b ^ I and found to be well borne out. 

A corresponding analysis for Sabiq) (not shown) finds 
the same agreement between the PRISM theory and the 
simulation data. These results indicate that a complete 
description of the static structure of our polymer melt 
can be achieved using only the site-independent inter- 
chain direct correlation function c(g) and the matrix of 
the single-chain structure factors Wabiq)- 
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FIG. 8: Comparison of the site-site inter-molecular pair- 
correlation functions hab{q) determined from the simulation 
data at T = 0.47 (circles) and from the PRISM theory (grey 
solid lines) [i^l . Some curves are shifted vertically for clarity. 
Note that h2i{q) = /12 10 (g) in the PRISM theory (see text for 
details) . 



C. Static Structure involving the Center of Mass 

1. Structure Factor of the Center of Mass 

At low temperature, the motion of the CM slows 
down similarly to that of the monomers [23|- For the 
monomers, we interpreted this behavior as a consequence 
of the tighter packing in the nearest-neighbor shells. If 
a tighter packing of the polymer coils and concomitant 
"coil-caging" was responsible for the sluggish dynamics 
of the chains, one could expect to find the signature of 
a stronger packing in the CM-CM structure factor Sc- 
Alternatively, polymer MCT |T3, Hlj suggests that the 
slowing down of the CM originates from the segment dy- 
namics which "enslaves" the CM motion. For this view to 
apply, little or no variation with T of the CM correlations 
is required, as the inter-molecular segment correlations, 
Pnih{q) (Fig. drive the (segment) caging. 

Recently, the spatial correlations between the CM's 
have been addressed in several studies [I^lisLli^lsollsil l. 




q 

FIG. 9: Static structure factor of the CM, Sciq), at all sim- 
ulated temperatures, i.e., T = 0.46, . . . , 1 (solid lines). The 
dashed line shows Eq. 1371 calculated from the simulation 
data for u;m,c(9), w{q) and h{q) at T = 1. The dotted hor- 
izontal line indicates the limit of Sciq) for g — > [Eq. 138ll 1. 
kBTpmKT was read off from the small-g behavior of S{q) 
(see Fig. The dotted vertical line shows 1/Rg — 0.69 
(iig ~ 1.45). 

References 0; Oil point out that deviations of the CM 
motion from free diffusion, observed for displacements 
smaller than the chain size already at high T, could be 
caused by inter-molecular interactions between the cen- 
ters of mass. Approximately, these interactions are given 
by the potential of mean force [s^, i.e., by the CM-CM 
pair distribution function. 

This function was also discussed in the context of de- 
vising efficient coarse-grained simulation models for semi- 
dilute polymer solutions and polymer-colloid mixtures 
,49 ■ 5Q]. In these studies, the attempt is made to repre- 
sent the polymer coils as soft, penetrating spheres. The 
spheres interact via an effective pair potential derived 
from the CM-CM distribution function. For this distribu- 
tion function Ref. [sj suggests a PRISM-approximation 
which relates Sc{q) to the structure factors of the 
monomers. In Appendix ^ we sketch the main ideas 
of this approach and discuss the validity of the under- 
lying assumptions for our polymer melt. Here, we only 
compare the result of the calculation, 

Sc{q)^l + ^ '''^ff pM, (37) 
N wyqY 

with our simulation data [w^,c{q) is defined in Eq. (|20|) ]. 

Figure|51shows Sc{q) for all investigated temperatures, 
together with Eq. 1)3 7|l. Starting from a small value {kt = 
isothermal compressibility) , 

(n^) — (n)^ 1 
lim Sciq) = ]—\ = -T^f^BTpmUT , (38) 

the structure factor increases toward a small peak before 
it approaches the ideal gas value 1 without any further 
oscillations. Although this peak is indicative of some 
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FIG. 10: Wm,c{qf /w{qf versus q at T = 0.46 (solid line) 
and T = 1 (dashed line). As w{q) does not depend on temper- 
ature (see Fig. the figure indicates that the monomer-CM 
correlation Wm,c is also independent of T. The dotted line 
shows the Gaussian approximation Eq. 13911 . The dotted ver- 
tical line indicates q = 1/-Rg — 0.69 (-Rg ~ 1.45). 



preferred distance between the chains, the effect is very 
weak. Furthermore, Sciq) is independent of T. This 
imphes that the sluggish dynamics of the CM is not re- 
lated to a tighter packing of the chains at low T. Rather 
it should be interpreted as a consequence of the slow- 
ing down of the monomer motion, which, due to chain 
connectivity, entails the glassy behavior of the CM. 

As observed in Ref. 51], we also find that the PRISM- 
approximation l(T7|) provides a good description of the 
simulation data, except for q close to the peak position, 
where Sc{q) is slightly underestimated. Thus, we may 
use Eq. (|37|l to interpret the finding that Sq (<z) does not 
change on cooling. Equation (|37|l contains the prefactor 
w-m,c{<lf' / w{qY , which depends on intra-chain correla- 
tions only, and is thus independent of T for our model. 
Figure [TUI compares Wrn.c.iQ)"^ IwjqY with the Gaussian 
approximation [see |53| and Eq. (PSl) ] 



w 



m,C 



(9) 



WYi{q) 



2(e 



- 1) 



(39) 



Quantitatively, Eq. (|39|l is not very accurate, presumably 
because our chains are too short. (The results for semi- 
dilute solutions of long chains obtained in Ref. [sj appear 
to agree better with Eq. H39|l .) Qualitatively however, 
the Gaussian approximation reproduces the simulation 
results. It starts at 1, has a maximum around q ~ 1.5, 
and vanishes for g > 5. Thus, the factor ■w„i,c{q)'^ /w{q)'^ 
eliminates the contributions coming from the local liquid- 
like structure of the melt (i.e., from h{q), see Fig.0} and, 
along with that, a possible dependence of Sciq) on T in 
our model. 



2. Correlation between the monomers and the CM 

Figure ITTI shows the site- resolved monomer-CM struc- 
ture factors. We see that the intra-chain contribution 



FIG. 11: Site-resolved structure factors resulting from 
monomer-CM correlations [Eqs. 117H19t 1: Wa,c{q) (intra- 
chain part), pha,c{<l) (inter-chain part), and Sa,c{<l) (^'U 
chains). These structure factors are (almost) independent 
of T. The data shown were obtained at T = 0.47. Circles 
indicate the Gaussian approximation [Eq. 14011 1 for Wa,c{q) 
at a = 1 (chain end) and a = 5 (middle monomer). The 
thickness of the lines and the symbols increases from a — 1 
to o = 5. 



Wa.cil) decays more slowly for the middle monomer 
(a = 5) than for the end monomer {a — 1). This observa- 
tion may be rationalized by the Gaussian approximation 

MM 



w 



a,C 



(q) = exp ■ 



1-3- 



N 



(40) 



Equation (|40|) is symmetric under a N — a. That 
is, chain ends are indistinguishable. The argument of 
the exponential is a parabola with a minimum at N/2. 
Thus, in qualitative agreement with the simulation data, 
^aci^) decreases more slowly with increasing q for the 
middle monomer than for the end monomers. In real 
space, this implies that the middle monomer is on average 
closer to the CM than the chain end ^5^. 

Figure 1111 also shows the pair-correlation function 
hafi{q) and the sum of intra- and inter-chain contribu- 
tion, Sa^cil)- Sa,c{<l) is related to the probability of 
finding a monomer at a certain distance from the CM 
of some chain, while ha^cil) measures this probability 
if the CM belongs to a different chain than that of the 
monomer a. Qualitatively, pha^ciq) appears to be the 
mirror-image of Wa^cil) with respect to the q-axis so that 
Sa.ciQ) is small. This agrees with the naive expectation 
that there is little correlation between the positions of 
the monomers and the CM's. However, Sa.cil) is not 
completely structureless. It exhibits a maximum for the 
middle monomer (a — 5), but a minimum for the end 
monomer (a = 1). Quite surprisingly, we find a posi- 
tive correlation of the middle segments and the CM's. 
On average, the probability of finding a middle (an end) 
monomer around the CM of a chain is increased (de- 
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FIG. 12: Comparison of the monomer-CM inter-molecular 
pair-correlation function ha,c{q) determined from the simu- 
lation data at T = 0.47 for a — 1 and a — 5 (circles) and 
from the PRISM theory (Eq. solid lines). The data for 

hi,c{q) are shifted vertically for clarity. 



creased) relative to random packing. Intra-molecular cor- 
relations are thus canceled by inter-molecular ones only 
at large distances (small q). At intermediate distances 
the intra-chain density distribution is either too little or 
too strongly compensated by the surrounding polymers. 
By averaging over all monomers along the backbone of 
the chain this site-dependence of Sa,c(,<l) is suppressed 
to a large extent. This means that a PRISM- like theory 
using the monomer averaged 5'm,c(9) (see Fig. I16f) only, 
could underestimate the monomer-CM coupling. 

However, this does not imply that the PRISM the- 
ory cannot be applied to explain the site-dependence of 
ha,c{q)- Equation H37|) results from the assumption that 
the CM may be treated as an additional, non-interacting 
site in the PRISM approach. That is, the monomer- 
CM- and the CM-CM direct correlation functions are 
supposed to vanish; only c{q) is kept. Using this as- 
sumption and 5*111, c/S' — "Wm.c/w (see Appendix IXjl we 
find from Eq. if^ 



Wm,c(g) 



r N 



^Waxiq) 



L x = l 



h{q) 



(41) 



Figure [T^ illustrates that Eq. H41|l is in good agreement 
with the simulation data. This allows two conclusions: 
First, Fig. 1 121 emphasizes again that the structural prop- 
erties of our model may be understood in terms of the 
site-independent inter-chain direct correlation function 
and site-dependent intra-chain structure factors. Second, 
the finding that ha,ciq) does not depend on temperature 
for our model (Fig. Ill|) may be explained by the same ar- 
gument put forward for Sciq)- It is related to the intra- 
chain contribution Wai^c/w{q) which is independent of T 
and suppresses the temperature dependence of h{q) for 
wave vectors around (/max (see Fig. I1U|I . 
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FIG. 13: Comparison of the three-monomer structure factor 
(thin lines) for the melt, Si{q,q,q), and for the polymers, 
'W3{q,q,q), with the respective convolution approximations 
(thick fines), ^(g)^ and w(g)^ at T = 0.47. The simulation 
results for the triple correlations are not smoothed. The lower 
statistical accuracy of 5*3(5, l) compared to S{q)'^ is clearly 
visible, especially at large q. 



V. RESULTS ON THREE-PARTICLE 
CORRELATION FUNCTIONS 

Recently, triple correlations in simple and network 
glass- forming liquids have been investigated [23|. This 
study shows that, while the convolution approximation 
(Eq. H24|) with C3 = 0) is very good for simple liquids, 
it fails to provide an accurate description of the cage 
structure in silica. As silica is a network-forming liq- 
uid, nearest-neighbor bonds make an important contri- 
bution to the local structure in the liquid. This is simi- 
lar to the chain connectivity in a polymer melt. So, we 
investigate the importance of triple correlations for our 
model by comparing the three-monomer structure factor 
with its convolution approximation for selected subsets 
of {q,k,p). 

Figure 1131 presents the three-monomer structure fac- 
tor of the melt, S^^q, q, q), and of the chains, ^3(9, q, q), 
for the choice that the three vectors q, k, and p form 
an equilateral triangle characterized by the length of its 
side q. We find that the convolution approximation pro- 
vides a very good description for the triple correlation 
of the polymers, the amplitude of the oscillations being 
slightly underestimated, however. For q < 20, S3{q^q,q) 
is equally well represented by the convolution approxi- 
mation, expect for a sharp dip at g « 6.3, revealing some 
anti-correlation at this g-value. For g > 20 the interpre- 
tation of the data is diflflcult due to the high noise level. 
In this region, we find that S3{q,q,q) is systematically 
larger than S{q)^ and even stays above unity, the the- 
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FIG. 14: w;3(i7, g,p = i7\/2(l — cos(/3)) (top) and SsCg, g,p = 
51^2(1 — cos ip)) (bottom) versus cos for some selected q- 
values at T = 0.47. The simulation data for wa and Sz 
are represented by thin lines, the convolution approximation 
[Eq. II24II with C3 — 0] by thick lines. Note that the data for 
Szisi) at g = 7.1 are rescaled by a factor of 0.1. 



oretical large q-limit of both quantities. This difference 
must be attributed to insufficient statistics ^3_7] . Despite 
this proviso, the convolution approximation represents a 
fairly good description of the three-monomer correlations 
for the choice of wave vectors (g, g, q). 

In order to investigate the angular dependence of 
the triple correlations we follow a suggestion made in 
Ref. [35]. We determine 6*3 and W3 for the triple of mod- 
uli (g, fc = (7,p = q^2{\ — cos (/s)), i.e., for isosceles tri- 
angles with two sides of length q enclosing an angle if> 
[Eq. (|23l]. 

Figure ^1 shows the simulation results and the con- 
volution approximation as a function of cos (/? for vari- 
ous q corresponding to maxima and minima positions of 
the three point structure factors (cf. Fig ll3ll . As found 
before, the agreement between W3 and the convolution 
approximation is very good, except at g = 24.9 where 
the approximation yields oscillations that are absent in 
W3. Similarly for most q- values, ^3 and its convolution 
approximation are fairly close to one another. Barring 
q = 24.9, for which the quality of the comparison is 
hard to judge due to the noise in S'3, noticeable devi- 
ations are obtained for wave vectors close to gmax- This 
might suggest that the cage structure in the cold melt 
imposes triple correlations which are different than those 
predicted by the convolution approximation. To test this 
conjecture the statistics of the data should be improved 



considerably, which is currently hard to achieve [31 



VI. SUMMARY AND CONCLUSIONS 

We explored static properties of a supercooled, non- 
entangled polymer melt consisting of flexible chains. The 
temperatures studied range from the high-T, normal liq- 
uid state of the melt to the supercooled state close to, but 
above the critical temperature of MCT (Tc > Tg). Our 
analysis utilizes various structure factors characterizing 
spatial correlations, on different length scales, between 
the monomers, between the monomers and the CM's, 
and between the CM's. The main findings of our work 
may be summarized as follows: 

Due to the flexiblity of the chains in our model (e.g., 
no bond-angle or torsional potential), their conforma- 
tional properties remain essentially unchanged on cool- 
ing. Thus, all intra-chain structure factors, w{q), Wab{q) 
and Wa,ci<l)i only depend very weakly on T (Figs. lM)lll|l . 

However, not only the intra-chain, but also the inter- 
chain correlation functions involving the CM, Sciq) and 
ha,c{<l), are (almost) independent of T (Figs. I^ITT|l . We 
explain this finding by PRISM theory which relates Sq (q) 
and hafiil) to the intra-chain structure by a term con- 
taining Wm,c/w{q). This ratio is (almost) independent of 
T (Fig. I1U|) and suppresses the temperature dependence 
of the inter-chain correlations on the local scale of the 
nearest-neighbor shells. 

The CM-CM structure factor Sc{q) is fairly feature- 
less (Fig. O. Outside the small-g regime refiecting the 
low compressibility of the melt, Sciq) quickly approaches 
the ideal gas behavior, Sc{q) = 1- This shows that the 
chains are soft, interpenetrating objects whose spatial ar- 
rangement does not change on cooling. For our model, 
the slowing down of the CM motion, observed in previous 
works 24] , is thus not related to a tighter packing of the 
polymers at low T. 

Correlations between the monomers are reflected by 
the collective structure factor S{q). S{q) exhibits hquid- 
like oscillations which become more pronounced with de- 
creasing temperature (Fig.j^J. This signature of a tighter 
packing in the nearest-neighbor shells results from the 
inter-chain contribution h{q) (Fig.^. Thus, the sluggish 
dynamics of our model is driven by the nearest neigh- 
bors that are not directly bonded to each other. In this 
respect, our polymer melt corresponds to a simple glass- 
forming liquid. 

Another result supports this correspondence. The 
convolution approximation which factorizes the three- 
particle structure factors 5*3(9, fc,p) in the product 
S{q)S{k)S{p) is generally invoked in the mode-coupling 
theory for the glass transition. Several studies show that 
this a ppr oximation is well justified for simple glass for- 
mers [23, 13 , but not for structurally more complicated 
ones, such as ort/io-terphenyl or silica In our 

model we find, analogously to simple liquids, that the 
convolution approximation works quite well, except for 
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Ssiq, k,p) at intermediate q (close to (/max) and for cer- 
tain angles between the wave vectors (Fig. I14|l . To what 
extent these deviations could be important in a mode- 
coupling calculation is hard to estimate quantitatively 
(due the lack of sufficient statistics 

The analysis of the monomer-resolved structure factors 
shows that the intra-chain contribution Wabiq) depends, 
to a very good approximation, only on the distance |a — fo| 
between the monomers (a, b) along the backbone. Thus, 
chain end effects are not very important for the intra- 
molecular structure of our model. On the other hand, the 
inter-chain structure factor hab{q) depends explicitly on 
the monomer pair (a, 6). The site-dependence of hatiq) 
and of hafiiq) may be explained by PRISM theory which 
assumes that the direct correlation function is indepen- 
dent of the monomer index (Figs. I8I12I) . By calculat- 
ing the site-site and site-averaged direct correlation func- 
tions we test this assumption and verify that it repre- 
sents a good approximation (Fig.lsj. This shows that the 
structural properties of our model, even subtle monomer- 
monomer- and monomer-CM correlations, may be calcu- 
lated from the site-averaged inter-chain direct correlation 
function and the site-dependent intra-chain structure fac- 
tors, both of which are determined in the simulation. A 
similar agreement between PRISM theory and computer 
simulations of coarse-grained polymer models was also 
found in other studies jiM Isil 155. .56] . In a forthcoming 
article, we will exploit the results of the present work to 
propose and test a mode-coupling theory for the dynam- 
ics of our supercooled polymer melt. 
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APPENDIX A: PRISM APPROXIMATION FOR 
THE STRUCTURE FACTOR OF THE CHAIN'S 
CENTER OF MASS 

The recent work by Krakoviack et al. [HJ] uses PRISM 
theory to calculate the CM structure factor Sc{q) from 
monomer correlation functions. This approach was found 
to compare well with simulation data of (long) chains in 



semi-dilute solution. For polymer melts Krakoviack et 
al. did not test the theory and even mention the caveat 
that the correlation between the CM and the monomers 
could be different. Here, we suggest that their theory 
also yields a reasonable approximation of Sq (q) for poly- 
mer melts. In the following, we sketch the main ideas of 
the approach of Ref. f5ll| and test the basic assumptions 
against our simulation results. 

The starting point of Ref. [sj consists in introduc- 
ing the center of mass as an additional, non-interacting 
site. This step generalizes the functions h, c, w, and S 
[Eqs. H13l - I15l33|l ] to 2 x 2 matrices 

(Al) 

c= I r^-' I , (A2) 

(A3) 

S=( .^^'"^ %^-'CM^| ^^^^ 

Here, the indices "m,m", "m,CM", and "CM, CM" denote 
the monomer-monomer, monomer-CM, and CM-CM cor- 
relations, respectively. By definition, all matrices are 
symmetric. With the notation of Sect. IIII Bl we have the 
following identities 

hm,m — h , ft.m,CM — ^m,C , /lCM,CM — he , 
Win,m = W , Wm,CM = W^,C , WCM,CM = 1 , (A5) 
S-m,m — S , 5'm,CM = 'S'm.C , <S'cM,CM = Sc ■ 

The matrices of correlation functions are related to 
each other by Eqs. I|30I35|I . i.e.. 









/iCM,CM 




Cm, CM 




CCM,CM 


, m 






WCM,CM 








S'cM.CM 



pc(g)=w-i(g)-S-i(g), 
h(g) = w{q)c{q) [w(q) + ph{q)] 



(A6) 
(A7) 



Using Eq. HA5|) we find from Eq. (|A6|) 

1 1 

/^mCni,m 
PmCm.CM 
PCCM,CM 



+ 



1 



/N ScS^S^jN' 
1 



(A8) 



Note that pm = Np [Eq.Q]. In Ref. ^] it is assumed 
that 



CCM,CM — CcM,m — Cm, CM — , 



(A9) 



and only Cm.m 7^ is retained. This implies that the 
centers of mass interact neither with each other nor with 
the monomers. 
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FIG. 15: Direct correlation functions versus q at T — 0.47. 
The correlation functions Cm,m(<7), Cm,cM('7), and ccM,CM(q) 
are calculated from Eq. (IA8II . The results of this calculation 
are numerically reliable for q > 1 IE3 . Additionally, the figure 
shows the monomer-monomer direct correlation c{q) given by 
Eq. and the CM-CM direct correlation function cc{q) 

defined from Sc{q) by pcc(g) = 1 — 1/5*0 (g). The monomer- 
monomer direct correlation functions Cm,m(<z) and c{q) are 
indistinguishable from one another (for g > 1 H^]). For q — > 
0, pmc(g) and pcc(g) tend to the limits [see Eqs. 13313811 1: 
Npn,c{q ^ 0) = 1 - N/ksTpn^HT = pcc{q 0) -840 for 
T = 0.47). 



To test the validity of this assumption we calculated 
the right-hand side of Eq. (|A8I) from the simulation 
data and checked whether the direct correlation func- 
tions involving the CM are small compared to Cm^ni- Fig- 
ure 1151 shows that the monomer-monomer direct corre- 
lation function is indeed (much) larger than the other 
components of c, at least for g > 1 [53 ■ Thus, Eq. (|A9|I 
appears to be a reasonable approximation for our model 
in the g-range, where Sc{q) substantially deviates from 
its small-g limit, Eq. H38|) (for a more precise test see 
Fig.Hni). 

Inserting Eq. (|A9|) into the Ornstein-Zernike equa- 
tion for h [Eq. IjATp ] and using Eqs. HA5I13|I we re- 
cover Eq. H33() for c^.miq) (= c{q)). Furthermore, if we 
insert Eq. HA9|I into Eq. HA8|I . we obtain the relation 
S'm,c — Wni,cS/w (for a test see Fig. lT^ from the second 
and the third line of Eq. IjASp . Using this relation and 
Eq. lO in the first Unc of Eq. we find 



1 



1 W, 



m,C 



(AlO) 



Pmh{q) , 



N w{q)^ 

p. Kg r ■' 

tained in Ref. [U (see Eq. (16) of Ref. 



i.e., Eq. (|^ . Eq uation (|A10|) is the central result ob- 
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